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ABSTRACT 



The relationship between the bending moment and the 
through-plane shear force in the vicinity of a mechanical fas- 
tener at failure was determined. Experiments were conducted 
on 4-inch wide flat plate aluminum and graphite-epoxy composite 
speciments that modeled portions of a wing skin along a spar 
and along a rib. The composite specimens were either 8-ply 
or 16-ply balanced layups and were simply supported at two 
opposing edges and free along the other two edges. The fas- 
teners were pulled normal to the plates, and the maximum force 
at failure was measured for specimen lengths varying from two 
to six inches between supports. 

The aluminum plates failed by formation of a plastic hinge 
across their width and showed little sensitivity to through- 
plane shear. The 8-ply spar specimens cracked across their 
width and also were relatively insensitive to through-plane 
shear. However, failures of the rib specimens were confined 
to a region near the fastener, where the fastener pulled 
through the plate, and showed much greater sensitivity to 
through-plane shear. 

Two analyses were made; one for small elastic deflections 
of a thin orthotropic plate, and another for a beam in the 
elastic range. A mesh generator for a finite element model of 
the plate around the fastener was also developed for the 
computer program ADINA. 
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INTRODUCTION 



I . 



A . BACKGROUND 

In recent years considerable effort has been devoted to 
the development of new materials for the purpose of improving 
aircraft performance. One family of these materials is fiber 
reinforced composites. Many studies have shown that use of 
advanced composite materials in aircraft structures can result 
in significant weight savings due to their low density, high 
modulus character. Furthermore, the anisotropic nature of 
fiber reinforced composites challenges the aircraft designer 
to exploit their directional nature in order to realize even 
greater savings. 

Graphite-epoxy composite materials have recently been speci- 
fied for major structural elements in two U.S. Navy and Marine 
aircraft currently in advanced development stages. Specifi- 
cally, the wing skin, leading edge extension, trailing edge 
flap, rudder, horizontal stabilizer, and vertical fin are to 
be constructed of graphite-epoxy in the F-13, [1] . The under- 

lying wing substructure is to be constructed of metal in a 
traditional spar-rib arrangement. The second aircraft, the 
AV-8B, is the advanced version of the AV-8 Harrier. This 
aircraft will have the entire wing box, consisting of the skin 
and the torque box substructure, constructed of graphite- 
epoxy, [2] 
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Mechanical fasteners are specified to attach the wing skins 
to the underlying wing substructure for both aircraft to allow 
removal of the skins, [3] . 

Extensive tests under various conditions have been con- 
ducted to assess the risks associated with graphite-epoxy 
materials. For example, strength degradation due to moisture 
absorption, high temperature environments, galvanic corrosion, 
and fatigue have been determined for the service environment 
of both aircraft, [3] . These tests have shown that design 
strain levels for both the Harrier and the F-18 wing designs 
provide adequate safety margins for the detrimental effects 
of elevated temperature and moisture absorption to be encoun- 
tered in service. 

In addition to the ability of the composite structure to 
withstand design loads under a wide range of service conditions, 
some attention has been given to the vulnerability of air- 
craft composite structures in a combat environment. Tests 
have been conducted to study the response of advanced composites 
to ballistic impact damage and the concomitant reduction in 
strength. A discussion of these tests is given in Reference [4]. 

Much data are available and a great deal is known about 
the behavior of metals at failure. Theories of failure have 
been tested time and again and failure models have proven va- 
luable in predicting stress levels that lead to failure. With 
the advent of composite materials, attempts to model the behav- 
ior of composites and the failure mechanisms have led to some 
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understanding of failure modes. However, significant gaps 
exist in the literature. Of particular interest here is the 
failure of a composite structure at a metal fastener in a 
particular mode called fastener pull-through. 

B. FASTENER PULL-THROUGH 

Consider a composite plate or skin fastened to the under- 
lying structure with mechanical fasteners, as in the F-18 and 
Harrier wing designs. Primary considerations in selecting skin 
thicknesses include in-plane shear, tensile strength and criti- 
cal buckling stresses. Joints are primarily designed to trans- 
fer in-plane tensile and compressive loads in the skin to the 
underlying spar-stiffener-rib structure. 

Three different failure modes at a joint are shown in Figure 
1. These are bearing, shear -out, and tensile failure. These 
failure modes are discussed in detail in Reference [5] . A 
fourth type of joint failure is illustrated in Figure 2. In 
this failure mode the fastener is pulled through the plate by 
a large through-plane shear force. Although very little atten- 
tion has been given to this type of failure in metals, this 
failure mechanism may be significant in fibrous composites 
under certain loading conditions . 

For example, through-plane shear stress exists in the 
skin of pressurized fuel tanks. Strain levels are usually 
small under normal design conditions, and the through-plane 
shear stresses at the joint are insignificant. However, if the 
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aircraft takes a hit in combat, circumstances may arise which 
subject joints around the fuel tanks to unusually high through- 
plane shear stresses- The extent to which the joint is able to 
carry these stresses may have significant impact on the vulner- 
ability of the aircraft. 

C. HYDRAULIC RAM 

Consider the problem of an integral wing fuel tank which is 
partially full of fluid when impacted by a high speed projectile. 
The kinetic energy of the projectile is transferred to the fluid 
in the tank as it passes through, [6] . The result of this en- 
ergy transfer has been characterized by three phases. The 
first phase is the propagation of a shock through the fluid. 

The second is the formation of a vapor filled cavity behind the 
projectile. The third phase is an oscillatory phase during 
which the cavity expands and contracts as the system returns 
to an equilibrium state. These phenomena, illustrated schema- 
tically in Figure 3, have been referred to as the hydraulic 
ram effect. 

The intense fluid pressure created by the hydraulic ram 
has a two-fold effect on the tank wall. First, the pressure, 
acting perpendicular to the wall, creates a moment distribu- 
tion through the plate. Secondly, the through-plane shear 
forces that arise must be reacted by the fasteners. This 
situation is shown in Figure 4. 
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Reference [7] describes hydraulic ram tests that were 
conducted on fluid containing aluminum tanks and one manifesta- 
tion of the failure of the tank walls was an "unzipping" of 
fasteners. Although the fasteners that were used in these 
tests were not typical aircraft fasteners, questions were 
raised about such connection failures. If aluminum plates 
and rivets failed in this through-plane mode, then it may be 
possible that composite plates will fail in the same way. 
Although composites are generally pound-for-pound stronger than 
aluminum in-plane, they may be considerably weaker in the 
through-plane direction. In addition to the bending moment 
and through-plane shear induced by the hydraulic ram loading, 
aircraft skins will be in tension, compression, and in-plane 
shear due to maneuvering when the hit occurs. This combina- 
tion of stresses can lead to a significantly lower failure load. 

D. OBJECTIVES 

Several aluminum and graphite-epoxy plates have been 
tested to determine the sensitivity of failure to the com- 
bination of through-plane shear and bending, moment due to 
hydraulic ram loading. Specifically, the objectives were to 
determine the force required to pull a fastener through the 
plates for a range of values of the bending moment at the 
fastener. A failure curve relating pull force and bending 
moment similar to a yield curve for limit analysis was sought. 
Since reinforcing fibers do not run in the through-plane 



direction in composite plates, it was expected that through- 
plane shear strength would be low compared to aluminum plates. 

Macroscopic failure mechanisms in both materials were 
observed. In order to verify experimental observations, two 
analyses were made. One analysis was for small, elastic de- 
flections to determine the magnitude of the errors in moment 
distribution associated with treating the specimens as beams 
instead of plates. Another model applies to beams in the 
plastic range. Additionally, a mesh generator for use in 
the elastic-plastic finite element code ADI NA was developed. 



16 



II. EXPERIMENTAL APPROACH AND SET-UP 



A . THE MODEL 

From Figure 4 it is apparent that the hydraulic ram 
loading on a plate creates a complex internal stress dis- 
tribution throughout the plate and particularly at the 
fastener. To further complicate matters, hydraulic ram is a 
dynamic problem involving large pressures applied over a 
short time interval. In order to facilitate the experimental 
and analytical approach, some simplifying assumptions were 
made. First, the dynamic effects of the hydraulic ram were 
not addressed, although it is known that they are of great 
importance in inducing catastrophic failures at fuel tank 
boundaries. Secondly, it was assumed that there would be no 
in-plane stresses present. Again, it was recognized that these 
will play an important role in joint failures. These assump- 
tions are justified here on the basis that this is a first 
study in an ongoing program. Future studies will consider 
these effects. 

Figures 5A and 5B illustrate the modeling of the skin 
with the test plates. The underlying structure may be envi- 
sioned as either a rib. Figure 5A, or a spar. Figure 5B. The 
difference has importance only with respect to the composite 
plates. For example, the zero degree ply direction was as- 
sumed to run parallel to the main wing spars, [1] . 
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Figure 5C shows the loading set-up to simulate the com- 
bination of bending moment and shearing force that exist at 
the fastener due to the hydraulic ram loading. The edges at 
y = ± b/2 were assumed to be free and those at x = 0 and 
x = a were assumed to be simply supported. The free condition 
was assumed because, near the rifcs and spars, moments due to 
the hydraulic ram loading would be smaller in the y direction 
compared to the moments in the x direction. 

B. TEST EQUIPMENT 

In order to load a specimen as suggested by Figure 5C, 
a devide was constructed in the Naval Postgraduate School 
Machine Shop. The simple support condition was provided by 
triangular knife edges fashioned from steel bar stock. A 
means was provided by which the supports could be positioned 
from one to five inches from the center of the specimen in 
one-inch increments. By moving the supports outward, the 
moment at the fastener becomes larger for a given pull force. 
Thus, supports close together provide a relatively larger 
pull force and smaller moment whereas supports farther apart 
give a relatively smaller pull force and larger moment. The 
supports were attached to a rigid steel base to facilitate 
attachment of the device to the moveable head of a Riehle 
200,000 lb. testing machine. The device is shown in Figure 6 
with an aluminum plate in place for a test pull. The fasteners 
were passed through a hole drilled in the specimens and at- 
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tached to a Clevis device which in turn was attached to a 
Baldwin SR-4 5000 lb. capacity load cell. 

The load cell was connected to a locally constructed 
Wheatstone bridge. A power supply and digital voltmeter 
were connected to the bridge circuit in the usual way and cali- 
brated to read an integer multiple of the load applied at 
the fastener. Bridge output was also monitored on a strip- 
chart recorder calibrated to read pounds of pull. The strip- 
chart recorder was most useful in providing a time history of 
the applied loading and assisted in recording significant 
events which occurred during the loading. Some aluminum 
specimens were fitted with strain gauges , which were wired to 
a Wheatstone bridge switching device. Gauge output was moni- 
tored by a second digital voltmeter calibrated to read micro- 
inches per inch. An additional channel of the strip-chart 
recorder was available to monitor the strain gauge output 
selected at the manual switching station. Figure 7 shows the 
entire experimental facility. 

C. TEST SPECIMENS 

1 . Aluminum Specimens 

The aluminum specimens were machined from sheet stock 
to a four- inch nominal width. Lengths of the specimens were 
varied from seven to fourteen inches, depending on the in- 
tended location of the supports. A hole was machined in the 
center of the specimen to provide a nominal two-thousandths 
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inch interference fit for the fastener as specified in 
Reference [8] . 

The fasteners that were used for the tests are the 
Hi-Tigue Hiloks that have been specified as skin fasteners 
for portions of the F-18. These HLT-328-8-18 fasteners 
are often specified for applications which require high fatigue 
resistance. Their unique Hi-Tigue feature, however, was not 
germane to the tests. Figure 8 shows the dimensions of the 
fastener . 

2 . Composite Specimens 

The composite plate specimens were fabricated by the 
author from Hercules AS-3501-6 twelve-inch pre-preg tape. 

The tape was a nominal ten mil thickness with low resin con- 
tent. Five 16-inch x 16-inch plates were constructed using 
either 8-ply [0/±45/90] s or 16-ply [0/±45/02/±45/90] s balanced 
layups and were prepared in accordance with Reference [9] . 

The composite preparation equipment is shown in Figure 9. 

Four-inch wide test specimens were cut from the 16- 
inch square plates after post-curing. Some 8-ply specimens 
were prepared such that the 0-degree ply was perpendicular to 
their axis (the spar specimens) . The remaining 8-ply and the 
16-ply specimens had the 0-degree ply parallel to the axis 
(the rib specimens) . The plates were drilled and fitted with 
the same fasteners as the aluminum specimens, although inter- 
ference fitting is not specified for this type of material. 
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Room temperature material properties from Reference [5] 
are tabulated along with computed flexural rigidities in 
Table I. 
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III. ANALYTICAL AND COMPUTER MODELS 



A. ANALYTICAL PLATE MODEL 



The test specimens are idealized as thin plates with 
two edges simply supported and two edges free, as shown in 
Figure 10. The deflections and moment distribution can be 
approximated by assuming the load to be uniformly distributed 
over the shaded rectangle of area uv . The governing equation 
for a balanced, symmetric laminated plate is given by [6] . 

D ll w 'xxxx +2(D 12 +2D 66 )w 'xxyy <!> 

+D z2 W 'yyyy =q<x ' y) 

where w is the plate deflection, q is the applied pressure 

load and subscripts denote partial derivatives. D , D , D 

11 12 22 

and D are stiffness coefficients. In the special case of 
6 6 

isotropy we have. 



D = D = D; D = uD; D = D(l-u)/2 (2) 

1 1 2 2 1 2 6 6 

The solution to Equation (1) for the portion of the plate 
PTSt shown in Figure 10 is given by [9] . 

w - ? { (a+A cosh (nury/a) + (3) 

+B m ( m '’ T -y/ a )sinh (m-n y/a ) } sin (m-rrx/a) 
where , 

a m = 4Pa^ (-1) / 2 sin (4) 

Ixx Z 3. 

5 5 _ 

t m D uv 

i l 

and a is the span of the plate. 
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Considering the unloaded portion of the plate beyond 
line ts , assume a deflection surface of the form. 



w 1 =E {A^coshfa Try/ a) + fin ny /a) sinh(irmy/a) 



m 



(5) 



+ C^.sinhfa 7ry/a ) + D ^ Try/a ) cosh fa.Try/a ) } s in(m nx/a ) 

The six constants A , B , A 1 , . ...D 1 , must satisfy the 

m m m m ■ L 

free edge boundary conditions at y = b/2, and continuity 
along line ts . Applying continuity conditions at y - b/2 
gives 



W = W 1 ; w , = W 1 , ; W , =h j , ; W , =W 1 , 

' x x' 'xx 'xx' 'xxx 'xxx 

It is shown in Reference [9] that 



A -A 
m 



m = a (y sinh 2y -cosh 2y ) 
m ' m ' m ’ m' 



( 6 ) 



(7) 



B -B 1 = a /2 (cosh 2y ) 

mm m' ' m 

C 1 = a (y cosh 2y -sinh 2y ) 
m m ' m ' m ' m y 



vl 



m 



= a m sinh 2 T m )/ 2 



where 



ym 



mirv 



m 



4a 



( 8 ) 



Two more equations are available from the boundary conditions 
at y=b/2 in order to solve for the six constants. Along the 
edge y=b/2 the moment M and the shearing force must 



vanish. Hence 
M 



- 


-(D w,„„ 


+ D w , „„ 


y=b/2 


22 yy 


12 xx 


1 


~(D w. 


+ (D +2 D 


y=b/2 


22 yyy 


1 2 



( 9 ) 



6 6 



xxy 



23 



Substitution of Equation (4) into Equation (9) gives, 



A* 

{ B . m > 

m 



4 Pa 4 (-1) O-D/ 2, 



m 5 ir 5 D uvK 
i i 



( 1 -p )a cosha -(p~-3)sinha (p -l)a sinha -2 cos 
^ 2 m m 2 m i m m 



ha 



m 



Jp 2 -l)sinha m 



Cl - P )cosha 



m 



/ \ 

-C 1 C 1 " P ) s inha -D 1 ( C 1 " P )a cosha +2 sinha 1 

m^ m m v ^ i ' mm 



nr 



-C 1 fl-p )cosha -D 1 ( fl-p la sinha -(p -3)cosha 1 

m^ 2 m m^ v 2 m m 2 nr 



where 



k = ( 1 -p 1 C 1 “ P )ot + [ 2 Cp ’ 1 ) + (p " 1 ) (p _ 3)] coshasinha^ 
v 1 v 2 m L 2 1 2 mm 



mTrb 

n - > P "" r\ 

1 r 22 



1 2 . _ (D 1 2 + 2 D 6 6 ) . = 

’ p 2 d 22 ’ p 3 



1 1 2 



m 2a ’ 1 D 22 ’ 2 D 22 ’ 3 D 11 

The quantity that is of particular interest here is the 
distribution of M along the line x = a/2. A plot of M /M 

A 

versus 2y/b is shown in Figure 11 for the aluminum specimen, 
and each of the three types of composite specimens for various 
values of a, where 
M = 

4 (ID 

B. PLASTIC LIMIT ANALYSIS FOR A BEAM 

If the aluminum specimen is considered as a beam,- then 
the total bending moment is given by Reference [12] as 

zcr(x) dz (12) 

t /2 



M = b 



J4 
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Where b is the width and t is the thickness of the beam. 



Making the usual assumptions that plane sections remain 

plane, the strain at any point is given by 

e = kz (13) 

Where k is the curvature of the middle surface, z = 0. For 

moments below a certain critical value, M , all stresses 

e 

are elastic and 

= Ee =Ekz (14) 

Substituting into Equation (12) gives 



M = ~~ b Ek t 3 



(15) 



If the moment is increased until at the outer fibers, 
z = ±t/2, is equal to the yield stress, a , the maximum 



(16) 



elastic moment is given by 

M = -i- b t 2 o 
e 6 y 

Since only the outermost fibers of the beam are yielded, the 
beam will continue to carry additional load through the cen- 
tral fibers until they too are yielded. The situation is 
shown schematically in Figure 12. The limiting moment is 
Case (D) of Figure 12 where the entire cross section is 
yielded and, 

1 



M 



° y bt‘ 



(17) 



o 4 

The moment in the beam is largest at the center of the 
span and its value is Pa/4. Substituting Equation (11) into 
Equation (17) , the maximum pull that the beam can carry, P , 



is given by 



p =a T fc2 
y y a 



(18) 



C. MESH GENERATOR FOR A FINITE ELEMENT MODEL 



The theory of matrix structural analysis is discussed at 
length in Reference [13] , and a discussion of various three- 
dimensional elements can be found in Reference [14] . The 
studies presented in Reference [14] by Clough demonstrate that 
hexahedral elements with nodes on the sides adequately repre- 
sent the bending of a simply supported plate. Clough used 
a single layer of 20 node elements and achieved the exact 
deflections. He showed that more complicated elements, such 
as curved tetrahedra, could perform as well, but that the 
formulation time was significantly greater. Further, it was 
shown that, while two-point Gauss quadrature rules gave exact 
results for rectangular prisms, three- or four-point rules 
were required for skewed elements in order to achieve the 
desired accuracy. 

In this study, 16 node isoparametric elements were chosen 
to model specimens. The objective of the finite element 
effort was primarily to develop a mesh generator which could 
later be used to analyze the laminated anisotropic plates. 

The finite element code used was AD IMA, developed by 
Bathe [15] in 1975. ADINA incorporates a library of four 
basic elements and allows selection of any combination of 15 
material models. The library material model that seemed 
best suited was the isotropic elastic-plastic model with strain 
hardening. Tensile tests were conducted on aluminum coupons 
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of three different alloys in order to accurately provide 
proper input data for ADINA . The results of the tensile 
tests were used throughout the investigation as baseline 
material properties, and the stress strain curves for three 
different aluminum alloys are shown in Figures 13A through 
13C. The bilinear approximation to the stress strain be- 
havior of the materials is indicated by the broken lines. 

Symmetry of the test specimens allowed some simplifi- 
cation of the model. The analysis includes one-quarter of the 
plate. The finite element model of the specimen is shown in 
Figure 14. The loads can be applied as concentrated loads at 
nodes along lines (A) and (B) of Figure 14. At x = 0, re- 
straints can be applied such that no displacement in the x 
direction is permitted. At y - 0, no displacement in the y 
direction is permitted. At x = a, the simple support condition 
dictates that no displacement in the z direction is permitted. 
The only other constraints involve the nodes at z = 0 and on 
the first ring of elements closest to the fastener. The assump- 
tion is made that the fastener is essentially rigid with respect 
to the aluminum. Consequently, no displacement in either the 
x or y directions is permitted for these nodes. 

A description of the mesh generator, and its use is in- 
cluded in Appendix A. Present limitations internal to the 
ADINA source code preclude the possibility of obtaining a 
successful nonlinear analysis due to array dimension. 
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IV. RESULTS AND DISCUSSION 



A. ALUMINUM SPECIMENS 

All of the aluminum specimens tested failed by forming a 
plastic hinge across the center of the span, as shown in 
Figure 15. None of the fasteners pulled through, and there was 
no visual indication that the through-plane shearing stresses 
had any effect on the failure mode. 

Table II contains the ultimate moment and pull force 
based upon the plastic limit analysis and the actual maximum 
moments, M^, and pull forces, P^, obtained during the experi- 
ments. Note that due to the fact that a = 4 inches, and 
Mf have the same numerical value. The ultimate moments ob- 
tained experimentally are seen to be somewhat higher than 
those predicted by the plastic limit analysis of section III. 
Part of the difference may be due to the fact that the theory 
applies to elastic-perf ectly plastic materials, whereas the 
materials tested (particularly the AL7075T6) exhibited some 
degree of strain hardening. Further, the limiting moments 
predicted by the theory are average moments, whereas the 
moment is not distributed uniformly across the beam, as shown 
in Figure 11. If the ultimate pull force is taken to be that 
force required to raise the moment M along the Line x = a/2 
to the ultimate moment, then the predicted values of pull, 

P f , more closely resemble those obtained experimentally. 

In any case, the tests indicate that these aluminum plates fail 
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at the joint in a plastic hinge mode. This hinge effect has 
previously been observed in hydraulic ram experiments, [7] . 

Figures 16 and 17 show a plot of the stress versus applied 
load for two of the locations which were provided with strain 
gauges on the 6061 aluminum specimen. The stress is normal- 
ized by the yield stress for the material, and the load is 
normalized by the theoretical P^ for the material and support 
conditions. The figures also contain the theoretical stress 
for a beam with a concentrated load of the same magnitude, 
and for the flat plate. For these tests, the simple supports 
were 4 inches apart. It can be seen that at the point x = 2, 
y = 1.040, the beam approximation does not predict stress levels 
accurately. The plate theory is considerably better. Further 
away from the centerline of the plate, at x = 2.735, y =.735, 
the plate theory overestimates the experimentally obtained 
stress levels. This suggests the possibility that formulation 
of the plate problem may not be accurate near the fastener. 

The assumption that the fastener head behaves like a uniform 
load acting over a rectangular area may not be refined enough 
to accurately predict the stresses in the vicinity of the 
fastener . 

B. COMPOSITE SPECIMENS 

The results of all tests on the composite specimens are 
presented in Table III. M u ]_t was determined f° r each of the 
five plates by applying a line load across the center of a 
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four-inch wide by seven-inch long specimen cut from each plate. 
The line load was applied by a steel cylinder attached to 
the Riehle testing machine, as shown in Figure 18A. The fas- 
tener passes through a hole in the cylinder and the hole in 
the plate and attaches to the 'dlevis. Thus, pulling the fas- 
tener causes the cylinder to bear against the plate along the 
axis of the cylinder. This eliminates the concentration of 
through-plane shear at the fastener. was determined 

for each 16-inch square plate by pulling a fastener through a 
two-inch by two-inch specimen restrained by a two-inch by six- 
inch by five -eighths - inch thick steel plate with a circular 
hole eighty- thousandths of an inch larger than the radius of 
the fastener head. The P ^ set-up is shown in Figure 18B. 

Two complete sets of experiments were run on the 8-ply rib 
specimens in order to establish that the results were not at- 
tributable to defects in the quality of the plates from which 
the specimens were cut. The mean difference was less than 7 
percent . 

Figure 19 is a plot of the experimental data in nondimen- 
sional form. The abscissa is the moment, M, normalized by 
M ^ . The ordinate is the pull force at failure, P£, normal- 
ized by P u £ . The extent or mode of failure is indicated in 
the figure. 

Figures 20A-20L consist of photographs of the 8-ply spar, 
8-ply rib, and 16-ply rib failed specimens. The sequence of 
photographs for each group of specimens is from small moment 



and mostly through-plane shear, to small through-plane shear 
and mostly moment loading conditions. The extent of the failure 
zone is indicative of the sensitivity of the specimen to through- 
plane shear. Note that the failure zone for specimens loaded 
with relatively small moment is primarily confined to a re- 
gion around the fastener hole. As the moment was increased 
(by lengthening the span) this failure zone extended in the 
y direction until it covered the entire width and resembled the 
failure of the line loaded specimen. 

Figure 21A is a plot of the data from a different point of 
view. The abscissa in this plot is the theoretical elastic 
local moment at the fastener head, , normalized by 

Figure 21B is a plot of the theoretical elastic moment at the 
edge of the plate at the center c£ the span, M e ^ CTe , normalized 
by . The theoretical values were obtained from Figures 

11A and 11B at y = 0.4 inches. In both cases the ordinate is 
the pull force at Failure, P^, normalized by . 
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V. CONCLUSIONS 



A. ALUMINUM SPECIMENS 

It was not possible to pull the HL-328-8 fasteners through 
the aluminum test plates. The failure of the aluminum joints 
was manifested entirely by the formation of a plastic hinge 
and the resultant inability of the specimen to carry further 
load . 

There is essentially no comparison between the way the 
aluminum specimens and the composite specimens failed when 
subjected to the combined effects of moment and through-plane 
shear. Failure of the aluminum plates was insensitive to the 
through-plane shear force induced by the head of the fastener. 

The assumption that the plates behaved like beams when sub- 
jected to a force on the fastener does not appear to be valid, 
with the exception of the special case where moment (or support 
spacing) is very large. It was shown from plate theory that 
the effect of the concentrated load upon the moment distribution 
is more significant where support spacings are small. The ap- 
proach taken here, therefore, requires some refinement before 
it can accurately predict stress levels. However, use of a 
finite element model should allow designers to more accurately 
predict stress levels that lead to failure. A wide variety 
of loading conditions and fastener geometries may be studied 
in this manner. 
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B. COMPOSITE SPECIMENS 



It is evident from Figure 20 that the weakly-oriented 
composite spar specimens, 52 and 53, were largely insensitive 
to through-plane shear. The location of the failure points on 
Figure 19 is also indicative of a lack of sensitivity since 
these specimens failed at a load that was very near the ulti- 
mate moment. However, the rib specimen showed a significant 
sensitivity to through-plane shear. The thickest specimens 
showed the greatest degree of sensitivity, and this result may 
be cause for concern. 

The manner in which pull-through failure occurs is complex. 
It is suspected that the ultimate moment is attained over an 
area near the fastener, causing failure of the specimen in 
that vicinity. However, the distribution of moments over the 
centerline is such that a portion of the plate outside this 
failure zone maintains sufficient strength to carry further 
load. The through-plane shear forces must still be carried by 
the composite around the fastener. Thus, as the load is in- 
creased further, the fastener pulls through the specimen. In 
those specimens which were loaded with large moment and small 
pull-through, the distribution of moments is flatter along the 
centerline, since the supports are further apart. That is, the 
effect of the concentrated load induced by the fastener head 
is less pronounced. Conseuqnetly , once local failure around 
the fastener occurs, only a small increase in pull force is 



required to increase the moment to its critical value along 
the entire width of the specimen. Thicker and wider speci- 
mens, in which the ultimate moment is larger, may show greater 
sensitivity to through-plane forces, and therefore be less 
likely to fail across their entire width. It would be in 
thick, stiff materials where "unzipping," if it is to occur, 
would be most likely. Current wing design concepts for por- 
tions of the F-18 wing call for thicknesses from .3 to about 
0.7 inches. Clearly, as thickness gets larger, the through- 
plane shear strength gets larger and there is a point where the 
actual through-plane shear loads can be carried by the com- 
posite . 
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VI. RECOMMENDATIONS 



A. ANALYTICAL METHODS 

A mesh generator has been developed to model the simply 
supported- simply supported- free - free plate loaded over a 
small region under the fastener head. This has been done 
using an elementary loading condition and an isotropic material. 
A great deal of understanding could be gained if the program 
were modified to include orthotropic elements. The current 
program allows generation of up to 16 layers of 16-node iso- 
tropic isoparametric elements. The node generation and ele- 
ment connectivity would remain unchanged. 

A further sophistication could be realized by 'treating the 
nodes near the fastener as fixed and applying pressure loads 
as predicted by the computer program Satans [16]. Satans has 
the ability to predict pressure distributions on a plate that 
arise in a fluid as a consequence of hydraulic ram from a 
specific threat projectile. Other modifications could include: 

1. The ability to change the geometry of the fastener head. 

2. The ability to take into account the flexibility of the 
fastener head. 

3. Inclusion of in-plane tensile, compression, and shear loads 
characteristic of maneuvering flight. 

B. EXPERIMENTAL PROCEDURES 

No data were obtained at the lower moment ranges because 
of the inability to move the support arms of the machine closer 



together than two inches. This design defect should be cor- 
rected by affixing the support arms to a base which permits a 
full range of motion, and does not restrict selection of dis- 
tances between arms to one inch increments. 

Investigation of the sensitivity of composites to a more 
complicated stress state which includes in-plane stresses is 
required. Further study is also required on factors which 
affect sensitivity, such as: 

1. Fastener head geometry variations, such as counter-sunk 
heads . 

2. Fastener-hole buffer materials and seals, such as 0-rings. 

3. Degree of interference fit or tolerance between fastener 
and hole. 

4. Effects of high cycle fatigue. 

5. Existing delaminations or other damage. 

6. Effects of softener strips. 



36 



APPENDIX A 



A. GENERAL 

The mesh generator is extremely simple to use. Where possi- 
ble, variable names are exactly the same as in Reference [15] . 

In addition to the source deck, only four data cards are re- 
quired to generate a data deck that can be used as input for 
ADINA. By simply removing the load cards, the deck can also be 
used with the preprocessor program PSAPl. The generator is 
specific and therefore limited in that it currently has the 
ability to generate nodes along 5 radial lines in one quadrant 
of a plate that has a hole in the center. Other nodes are less 
restricted. The group of elements closest to the fastener hole 
are fixed in size by fastener parameters that are input. The 
size of the rest of the elements is determined by the number of 
divisions desired along the line that lies at the center of the 
span of the plate. The capability currently exists to generate 
up to 16 layers of 16 node isotropic isoparametric elements. 
However, the source program ADINA requires minor internal modi- 
fication to accommodate the large volume of storage required by 
a non-linear analysis. A program is currently underway at the 
Naval Postgraduate School to provide ADINA with this increased 
capability . 

B. REQUIRED INPUT FOR MESH GENERATOR 

Input data for the generator requires only four cards. 
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1. Input Card 1 Format (1615) 



COLUMN 


VARIABLE 


1 - 


5 


IDOF(l) 


6 - 


10 


IDOF (2) 


11 - 


15 


IDOF (3) 


16 - 


20 


IDOF (4) 


21 - 


25 


IDOF (5) 


26 - 


30 


IDOF (6) 


31 - 


35 


NEGNL 


36 - 


40 


MODEX 


41 - 


45 ■ 


NSTE 


46 - 


50 


IPRI 


51 - 


55 


NPB 


56 - 


60 


LAYERS 


61 - 


65 


NEGL 


66 - 


70 


NLOAD 


71 - 


75 


NLCUR 


76 - 


80 


NPTM 


Input 


. Card 


2 FORMAT ( 2 0A4 


1 - 


80 


PTYPE 



DESCRIPTION 



Number of layers of elements 
desired . 



ALPHANUMERIC TITLE 

DESIRED ON 

ADINA AND PSAPI PRINTOUTS. 
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3. INPUT CARD 3 FORMAT (8F10.0) 



1 


- 


10 


PROP (1) 


Young's Modulus 


11 


- 


20 


PROP ( 2 ) 


Poisson's Ratio 


21 


- 


30 


PROP (3) 


Simple tension 


31 


- 


40 


PROP ( 4 ) 


Strain Hardening Modulus 
A Length of 


41 




50 




half span between simple 
supports 


51 


- 


60 


Rl 


Radius of Fastener 


61 


- 


70 


PMAX 


Maximum load desired 


71 




80 


PINC 


Desired load increment for 
load data cards. Different 
from DT which is load step 
internal to ADINA, and deals 
with the method of solution. 


INPUT CARD 


4 FORMAT (8F10. 


0) 


1 


- 


10 


R 


Radius of Fastener Head 


11 


- 


20 


T 


Thickness of specimen 


21 


- 


30 


DT 


Load step increment 
Allow to default (0.0) 


31 


- 


40 


T START 


Allow to default (0.0) 


41 


- 


50 


FAC 


Allow to default (0.0) 
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OOOOO CD OOOOOOo-g O Ui oooo 



G£NERATGR PROGRAM L I S TI NG** 3 ***#* 






DIMENSION ANGL(9) ,10(3,1500 ),X( 1500), Y( 1500 ), U 1500), 
*KEL( 3 2, 16),KG( 32 » 16) ,N0D( 16 ,3 2, 16) ,1 PN0D£( 2 ,3) , I DCF (6) 
*NFAR(20) , PROP (4),TIMV(20),RV(20),RV1(20),N0DE( 1500) 
READ( 2, 5) ( ID0F ( I ) , I =L , 6 ) , NEGNL , MCDEX , NST £ , IPRI,NPB, 

*L AY £RS , NEGL ,NL0A D , NLCUR , NPT M 
R£AD( 2*5) NDI V , I D C , ICON, IPRIC, MAX ES 
FORMAT ( 1615) 

R EAD( 2,6) PTYPE 
FORMAT ( 20A4) 

R EAD (2,7) (PROP! I ) , I = 1 , 4) ,A ,R 1 , P MAX, PI NC 
R EAD (2,7) R2,T, DT ,T START , FAC 
FORMAT ( 8F10.0) 



ASSIGN VALUES TO CONSTANTS AND SET DEFAULTS 



NPR=2*NDI V+3 
NPRM I D= (NPR + 1 )/ 2 
DELR=( A-R2) /FLOAT ( 2*NDI V ) 
N = 0 

NODE ( 1 ) =000 00 
NN=5*NPR+4*NPRMID 
NEL=4*( NDIV+1) 
NDIV1=N0IV+1 
PI =4. 0*ATAN( 1. 0) 

ANGL ( 1 )=0 .0 

ANGL ( 2 ) =ATAM .25 ) 

ANGL(3)=ATAN(.5) 

ANGL ( 4 )= AT AM .75) 

ANGL( 5) =PI/4.0 

ANGL ( 6)=ATAN( 4. 0/3.0) 

ANGL(7)=ATAN(2.0) 

ANGL ( 8) =ATAN( 4.0) 

ANGL (9 )=P 1/2.0 
LAYRS=LAYERS+1 
I PNODE ( 1, 1 ) =1 
I PNODE ( 2, 1 )=NN 
I F ( NPB. EQ. 2) GO TO 8 
I PNODE( 1,3) =NN + 1 
I PN0DE( 2,3 )=NN*LAYERS 
I PNODE (1,2) =NN*LAYERS+1 
I PNODE ( 2, 2 )=NN*LAYRS 
I F ( NLOAD. EQ.O ) NLOAD=23 
IF(NLCUR. EQ .0) NLCUR=2 
IF(NPTM.EQ.O) NPTM=NSTE 
M TYP = 1 
I ELD=0 
I ELX=0 
I ST = 0 
I 0 IRN= 1 
NPT S=NPTM 
N TF 1= 1 
NT F2 = 2 



GENERATE THE NODAL MESH, LAYER BY LAYER 
EACH ELEMENT WILL HAVE 16 NODES, 8 ON THE 
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C Y-Z UPPER AND 8 ON THE Y-Z LOWER SURFACE. 

C NO NODES ARE GENERATED THROUGH THE THICKNESS 

C OF THE ELEMENT. 

C 

DO 50 K=I , LAYRS 

XN=T * ( I .0-FL0AT( K- I) /FLOAT ( LAYERS) ) 

DO 40 J = 1 » 9 
THETA=ANGL( J) 

I F ( MOD ( J t 2 ) . EQ . 0 ) GO TO 21 
DO 31 1= 1 t NPR 
N=N+ 1 

IF(N.NE.l) NODE ( N)=N0DE(N-1 ) + l 
IF(N.EQ.l) NODE ( 1 ) =1 
ID ( 1 » N ) =0 
I 0 ( 2 » N ) =0 
I D( 3, N) =0 

IFIJ.EQ.l) I D( 3 » N ) = 1 
I F ( J.EQ.9) ID(2,N) =1 
I F( J.GE .5. AND. I .EQ.NPR) I D t I ,N) =1 
IF(K.NE .LAYRS) GO TO 15 
IFII.EQ.1) I D( 2 » N ) = 1 
IF(I.EQ.l)ID(3 t N)=l 

15 X ( N )=XN 
IF(I.NE.l) GO TO 16 
Y(N)=R1~CQS(THETA) 

Z(N)=R1*SIN(THET A) 

GO TO 30 

16 IFII.NE.2) GO TO 17 

Y ( N) = ( R1+(R2-R1)/2.0)*C0S(THETA) 

Z IN ) = (R1+(R2-R1) / 2.0) *SIN( THETA) 

GO TO 30 

17 I F ( I.NE.3) GO TO 18 

Y (N)=R2*C0S (THETA ) 

Z(N)=R2*SIN(THETA) 

GO TO 30 

18 I F( I. EQ.NPR) GO TO 19 

Y I N ) = I R2 + F LCAT ( I— 3)*DELR)*CQS(THETA) 

Z ( N ) = ( R 2+FLQAT ( I - 3 ) *D ELR ) * S IN ( THET A) 

GO TO 30 

19 I F( J.GT.5) GO TO 20 

Y ( N ) = A 

Z(N)=A*TAN(THETA) 

GO TO 30 

20 Y (N)=A*COTAN(THETA) 

Z ( N ) =A 

30 IF(J.EQ.9)Y(N)=0.0 

31 CONTINUE 
GO TO 40 

21 DO 33 1=1 »NPRMID 
N=N+1 

NDINC=2 

IFII. EQ.1 )NDINC=1 
NODE! N) =NODE (N-l ) +1 
ID! 1 » N ) =0 

I 0 ( 2 » N ) =0 
ID! 3,N) =0 

IFIJ. GE .5. AND. I.EQ.NPRMIDHDI 1»N)=1 
I F(K. EQ. LAYRS. AND. I.EQ.l ) I D ( 2»N ) = 1 
ID! 3,N) =ID( 2,N) 

X (N)=XN 

IFI I .NE .1 ) GO TO 22 
Y(N)=R1*C0S(THETA) 

Z (N) = R1 *S IN (THETA ) 

GG TO 32 

22 IFI I .NE . 2 ) GO TO 23 

Y I N ) = R2 *C OS (THETA) 

Z(N)=R2*SIN( THETA) 

GO TO 32 

23 IFI I .EQ.NPRMID) GO TO 24 

Y IN )= I R 2+FLOAT I I-2)*2.0*DELR)«C0S(THETA) 

Z ( N ) = I R2 + FL0AT I I - 2 ) *2 . 0*DELR ) * S IN I THETA ) 
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24 



25 

32 

33 
40 
50 
C 

C_. 

C 

C 

C 

C 

C 

C 

C 

C 

C_ 

C 



55 



56 



56 



59 



57 

99 

C 

C__ 

c 

c 

c 

c 

c 



GO TO 32 

IFIJ.GT.5 )G0 TO 25 
Y(N)=A 

ZIN)=A*TANI THETA) 

GO TO 32 

Y(N)=A*COTAN(THETA) 
Z (N)=A 

IF ( J.EQ.9)Y(N)=0.0 
CONTINUE 
CONTINUE 
CONTINUE 



GENERATE ELEMENT CONNECTIVITY. THIS IS 
DONE BY SETTING UP THE FIRST ELEMENT IN 
EACH LAYER, AND REFERENCING THE REST 
OF THE ELEMENTS IN THAT LAYER TO IT. THE 
TASK IS SIMPLIFIED TO SOME EXTENT BECAUSE 
ADINA HAS PROVISIONS FOR ELEMENT GENERATION. 



DC 99 1=1, LAYERS 
KE L ( 1, I)=( 1-1 )*1 6+1 
KG( 1, I )=0 

NODI I , i , I ) = I I-I ) *NN+1 

NODI 2,1,1 ) =NQD I 1 , l ,1 ) +2 

NOD I 3,1 ,1 )=N0D(2, 1,1) +NPR+NPRMID 

NODI 4,1 , 1 ) =NQD I 3 , 1 , I ) -2 

DO 55 L= 1 » 4 

NODI L + 4 ,1 , I )=NOD( L, 1, I )+NN 
CONTINUE 

NCDI9, 1 , I )=NOD( 1 , 1, I ) +1 

NODI 10,1,1) =NQD I 2 ,1,1 J+NPR-l 

NODI 11,1,1 ) =NQD ( 9,1,1 ) + NPR+ NPRM I D 

NODI 12, 1, I ) =N0 D ( 10, 1, I )-l 

DO 56 L = 1 , 4 

NOD I L + 12, 1 , I ) =NOD I L + 8 , 1,1 )+NN 

CONTINUE 

KOUNT= 1 

DO 57 J = 2 , N EL 

KCUNT=KOUNT + 1 

KGI J , I ) =0 

KELI J, I )=KEL( J-l, I ) + 1 

NODI 1 , J , I ) =NOD I 2 , J-l , I ) 

IFIKOUNT.EQ.INDIV+2) ) NODI 1,J,I ) =NOD( 4 , J-ND I VI , I ) 

NOD 1 2 , J , I ) =NOD I 1 , J , I ) +2 

NODI 3, J,I)=NOD( 2, J,I ) +NPR+NPRM I D 

NODI 4, J , I ) =NOD( 3 , J , I )-2 

DO 58 L=1 , 4 

NODI L+4, J, I ) =NOD I L,J, I)+NN 
CGNT INUE 

NODI 9 , J , I ) = NOD I 1 ,J,I)+1 
NODI 10, J, I ) =NOD I 10*J-1,I) + 1 

IFIKOUNT.EQ. (NDIV + 2) )NOD( 10 , J, I )=NOD( 10, J-ND I VI, I )+NPR 

NODI 11, J, I )=NOD( 3 , J, I )-l + N P RM T n 

NODI 12, J, I ) =NOD( 1 0 , J , I )- 1 +NPRMI D 

DO 59 L=1 ,4 

NODI 12+L, J, I ) =NOD I 8+ L , J , I ) + NN 
CCNT INUE 

IFIKOUNT.EQ.INDIV+2) )K0UNT=1 
MAXE L=KEL I J , I ) 

CONT INUE 
CONTINUE 



DIRECT OUTPUT TO THE APPROPRIATE DEVICE. FOR 
DATA CHECK FUNCTIONS THIS 'WOULD BE THE OFFLINE 
PRINTER. WHEN THE DATA IS SATISFACTORY, OUTPUT 
BE DIRECTED TO THE OFFLINE PUNCn. IT 1$ THEN RE 
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TO BE USED WITH THE ADINA JOB CONTROL DECK. TO 
USE IT WITH PSAP1 SIMPLY REMOVE THE LOAD CARDS 



C 
C 
C 
C 
C 

WRITE! 7 » 1 60 ) MA XE L » PTYPE 
160 FORMAT! IX, 15, 2X, • ELEMENT PLATE*, 15A4) 

WRITE!7,170)N, < I DOF! I), I = 1,6 1 , N EGL ,N EGNL , MODEX , NS TE , 
*DT,TSTART ,IPRI 

170 FORMAT! 15, 611, 14, 315, 2F 10. 0,15) 

WRITE! 7,180) 

180 FORMAT!///) 

WRITE ! 7 , 1 90 ) NPB, IDC 
190 FORMAT! 1615) 

WRITE(7,190)!( I PNODE ( I,J),I=I,2),J = i,NPB) 

DC 210 1=1 , N 

WRITE! 7,200) NODE <I),!ID<J,I),J=1,3),X<I),Y(I),Z!I) 

200 FORMAT! IX , I 4 , 1 X , I 4, 2 I 5 , 15X , 3 E 1 0 . 3 ) 

210 CONTINUE 

WRITE ( 7, 190)NLOAD,NLCUR,NPTM 
WRITE (7 ,190 ) ICON, IPRIC 
DO 220 1=1,20 
NPAR! I ) =0 
220 CONTINUE 
NPAR! 1 ) =3 
NPAR ( 2 ) =MAX EL 
NPAR! 3 ) =1 
NPAR! 7) =16 
NPAR! 10 )=3 
NPAR! 1 1 ) =3 
NPAR! 1 5 )= 8 
NPAR! 16 ) = 1 

WRITE! 7,230)(NPAR(I) ,1=1,20) 

230 FORMAT! 2014) 

WRITE(7,190) MTY P 

WRITE! 7 ,240) { PROP ! I ) , I = 1 ,4 ) 

240 FORMAT ( 8E10 .3 ) 

DO 270 1=1, LAYERS 
IPS = 0 

IFd.NE.l . OR. I. NE. LAYERS) IP S= I 
DO 260 J=1 , NEL 

WRITE! 7,25 0) KEL! J, I ) , I ELD , I ELX , I PS , MTY P , MA XE S, I ST, KG (J y l ) 
WRITE! 7,190) (NOD ( K, J , I ) ,K=1 ,8) 

WRITE! 7, 190)! NOD! K,J,I ) ,-K=9 ,16) 

250 FORMAT! 8 15, E10.3 ) 

260 CONTINUE 
270 CONTINUE 
C 

C 

C LOAD CARD GENERATION AND OUTPUT FOLLOWS 

C 

C _ 

C 

DO 280 1= 1 , NSTE 
T I MV! I ) =1-1 

R V 1 1 ) =P INC/ 80. 0* FLOAT (I ) 

IF(RV! I ).GT .(PMAX/80.0) )RV( I ) =P MAX/80.0 
RV1 ( I ) =RV ! I )/2 .0 
280 CONTINUE 

WRITE! 7, 190 )NTF 1 , NPTS 

WRITE! 7,240) (TIMV! I ) ,RV(I ), 1=1, NPTS) 

WRITE! 7,190)NTF2»NPTS 

WRITE(7 ,240)(TIMV ( I) ,RV1( I) , 1=1, NPTS) 

DO 300 1=1,9 
NCUR= 1 

IF! I.EQ.l .OR. I .EQ.9) NCUR = 2 
I F( MOD! I , 2) .EQ. 0) GO TO 290 
NLDl=LAYERS*NN+l+< 1-1 )*! NPR +NPRM I D ) / 2 
NLD2=NLD1 + 1 
NLD 3=NL D2+ 1 

WRITE! 7,3 10 )NLD1, IDIRN,NCUR, FAC 
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290 

300 

310 



WRITE! 7,31 0)NLD2 » IDIRN, NCUR , FAC 
WRITE ( 7 , 3 10 )NLD3, I DI RN , NCUR , FAC 
GO TO 300 
NLD 1=NLD1+NPR 
NLD2 = NL DI +1 

WRITE! 7 , 3 10 ) NLD 1 , I DI RN , NCUR , FAC 
WRITE! 7,310 )NLD2 » I DI RN , NC UR , FAC 
CONTINUE 

FORMAT! 3I5t F10.0) 

STOP 

END 
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E = 21.0 x 

l l 


10 6 psi 






E = 1.7 x 

2 2 


10 6 psi 






G = .65x 

1 2 


10 6 psi 






v = .21 

1 2 








v = .017 

2 1 








8-ply spar 


8-ply rib 


16-ply rib 


D 

l 1 


169.3 


634.4 


4815 


D 

1 2 


99.7 


89.7 


914.1 


D 

2 2 


634.4 


169.3 


1379 


D 

6 6 


112.1 


112.1 


1014 


D 

1 6 
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0 
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2 6 
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0 



High Strength Graphite-Epoxy material properties 



TABLE I 
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MATERIAL 





2024T3 


6061T6 


7075T6 


t (ins) 


0.25 


. 080 


. 050 


a (ins) 


4.0 


4.0 


4.0 


Theoretical 
Mo, (in-lb) 
Py , (lb. ) 


3520. 


272. 


175. 


Experimental 

M f,p f 

(in. lb.) (lb.) 


4800. 


380. 


246. 


Theoretical 
(plate) Mo, 
P 

y 


4928 . 


310. 


240. 



TABLE II 
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Figure I 




Pull-through failure of a 
mechanically fastened joint 

f 

Figure 2 
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Oscillation of cavity 



Figure 3 
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Pressure, stress and moment 
distribution induced by 
hydraulic ram 



Figure 4 
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test 



t es t 




(A) 








The model 
Fig u r e 5 
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Simple support device with 



aluminum specimen in place 



* 







Figure 6 
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Experimental facility 







HL 328-8-18 fastener used in the tests 
(Patent held by Hi-Shear corp.) 



Figure 8 
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Composite preparation equipment 





Test specimen idealization 



Figure 10 
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elastic 



elastic limit 



(A) 



CB) 



tfy 





plastic limit 



(C) 



(D) 



Stress distribution in a rectangular beam 



Figure 12 



T 



S 

N 



t 1 1 1 r 



□ 



s 

IS 

« 



□ 



in 

n 

h 

X 

P4 

ts 

W 

-1 

££ 



V 

P 

a 

□ 



□ 



is 

m 

w 

IS 

CM 

o 



x 







C3 

X 

• 

IS 



p 



□ 



s 

P4 

ts 



p 

J I I I I I L 



ca 


ts 


ts 


ts 


ca 


ts 


ca 


ts 


ts 


ts 


ts 


IS 


IS 


ca 


• 


• 


« 


• 


« 


• 


• 


r- 


in 


1/1 


X 


n 




— 



( ^.0 I * !Sd) p 



ts 

IS 

« 

J a 

ts 

B 

ts 



63 



Figure 13 A 
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Figure I3B 
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Figure I3C 
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Aluminum specimen 
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Figure 17 
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(A) a= 2", spar specimen 
P f =452 lbs. 



(B) 



a = 2", spar specimen, 

I ine load , P f = 5 30 lbs. 



Figure 2 0 
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1C) thr ou gh-plxine shear, (D) o=2”, 8-ply rib specimen, 

small moment, 8-ply, P = 605 Lbs. 

P f =113 5 lbs. 





(E) a = 3“, 8-ply rib specimen, 
P f = 582 lbs. 



(F) 



a=4", 8-ply rib specimen, 
P =410 l_b s. 



Fi g ure 2 0 (cont'd) 



7 3 










(G) a=5", 8-ply rib specimen, 
P f = 302 lbs. 



(H) a = 2", 8-ply rib specimen, 
line loaded, P f - 1235 lbs. 




(I) 



through-plane shear, 
small moment, 16-ply, 



P f = 2925 lbs. 



IJ) a = 3", 16-ply rib specimen, 
P = 1255 lbs. 



Figure 20 



7 4 



(cont’d) 













(K) a= 4", 16-ply rib specimen, ( L) a=5”, 16-phy rib specimen, 
P f = 1140 lbs. P f =1085 lbs. 




(M) a = 6”, I 6-ply 
P f =-985 lbs. 



rib speci men, 



( N) a= 4", 1 6-p I y ri b 
line loaded, P 



s p ec i m en, 
= 1780 lbs. 



Figure 20 (cont'd) 
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Figure 2 I A 
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Figure 2IB 
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